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Abstract 

Mott insulators with a half-filled band of electrons on the triangular lattice have been recently 
studied in a variety of organic compounds. All of these compounds undergo transitions to metal- 
lic/superconducting states under moderate hydrostatic pressure. We describe the Mott insulator 
using its hypothetical proximity to a Z2 spin liquid of bosonic spinons. This spin liquid has quan- 
tum phase transitions to descendant confining states with Neel or valence bond solid order, and 
the insulator can be on either side of one of these transitions. We present a theory of fermionic 
charged excitations in these states, and describe the route to metallic states with Fermi surfaces. 
We argue that an excitonic condensate can form near this insulator-metal transition, due to the 
formation of charge neutral pairs of charge +e and charge — e fermions. This condensate breaks 
the lattice space group symmetry, and we propose its onset as an explanation of a low temperature 
anomaly in k-(ET)2Cu2(CN)3. We also describe the separate BCS instability of the metallic states 
to the pairing of like-charge fermions and the onset of superconductivity. 



I. INTRODUCTION 



The organic superconductors have proved to be valuable systems for exploring strong cor- 
relation physics near the metal-insulator transition: it is possible to tune parameters across 
the metal-insulator transition by moderate pressure, while maintaining a commensurate den- 
sity of carriers. Of particular interest are the series of compounds which realize a half-filled 
band of electrons on a triangular lattice. In the insulator, three distinct ground states are re- 
alized in closely related compounds: (i) a Neel ordered state^in k-(ET) 2 Cu[N(CN) 2 ]C1 (ii) 
a recently discovered valence bond solid (VBS) stat d 2 * 3 * 4 1 in EtMe3P[Pd(dmit)2]2, and (Hi) an 
enigmatic spin liquid state™ 3 in fi>(ET)2Cu2(CN)3. Under pressure, all these compounds 
undergo a transition to a superconducting state. 

The nearest-neighbor 5 = 1/2 Heisenberg antiferromagnet on the triangular lattice is 
believed to have a small, but non-zero, average magnetic moment!- 8 * 9 * 10 * 11 ! However, pertur- 
bations such as ring and second-neighbor exchange are likely to be present near the metal- 
insulator transition, and these presumably yield the non-magnetic states noted above. This 
paper will examine the organic compounds from the perspective of a particular spin-liquid 
insulating ground state: the Z 2 spin liquid of bosonic spinons proposed in Refs. fT2lfT3l An 
advantage of this state is that it is a very natural starting point for building the entire phase 
diagram: theories of the transitions out of this Z 2 spin liquid state into the magnetically 
ordered (Neel) state^l and the valence bond solid stated have already been presented, thus 
realizing the three classes of insulators which are observed in these compounds. 

The core analysis of this paper concerns the insulator-to-metal transition on the triangular 
lattice^ out of the Z 2 spin liquid state. The spin-charge separation of the insulating spin 
liquid state (and the associated Z 2 topological order) will survive across the transition into 
the metallic state, and so this is a insulator-metal transition between two 'exotic' states. 
However, our theory generalizes easily to the other 'non-exotic' insulators with conventional 
order noted above: in these cases it is the conventional order which survives across the 
transition, rather than the topological order. Thus our study of the Z 2 spin liquid state 
offers a convenient context in which to study the insulator-metal transition, with a focus 
on the charged excitations in the simplest situation. Subsequently, we can include the spin- 
sector instabilities of the Z 2 spin liquid to conventionally ordered confining states; these do 
not modify the physics of the insulator-metal transition in a significant manner. 

Our perspective on the insulator-metal transition should be contrasted to that in dy- 
namical mean- field theory™ 3 (DMFT). DMFT be gins with a correlated metallic state, and 
presents a theory of the metal-to-insulator transition. Partly as a consequence of its focus 
on the metal, it does a rather poor job of capturing the superexchange interactions in the 
insulator, and the consequent resonating valence bond (RVB) correlations. Instead, we be- 
gin with a theory of RVB-like correlations in the insulator, and then develop a theory of 
the charged excitations; the closing of the gap to these charged excitations describes the 
insulator-to-metal transition. A theory of the closing of the charge gap in an RVB insulator 
has also been discussed recently by Hermele^, but for an insulator with gapless fermionic 
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spinons on the honeycomb lattice. 

For a half-filled band, our main result is that prior to reaching the metallic state, the 
insulator has an instability to a state with exciton condensatiot¥^, associated with the pair- 
ing of fermions of charge e with fermions of charge — e. As in the conventional theory^, 
the pairing arises because of the attractive Coulomb interaction between these oppositely 
charged fermions, while the charge gap to these fermionic excitations is still finite, but small. 
The excitonic insulator maintains the charge gap, but is associated with a breaking of the 
lattice space group symmetry which we shall describe. Eventually, Fermi surfaces of charge 
±e fermions appear, leading to a metallic state in which the excitonic instability (and the 
lattice symmetry breaking) can also survive. 

We propose that the above instability offers a resolution of a central puzzle in the prop- 
erties of the 'spin liquid' compound k-(ET) 2 Cu 2 (CN) 3 . At a temperature around 10K a 



feature is observecP^ in nearly all observables suggesting a phase transition into a mysteri- 
ous low temperature correlated state. An important characteristic of this compound is the 
indication from optical conductivity measurements^ that the charge gap of this insulator 
is quite small, suggesting its proximity to the metal-insulator transition even at ambient 
pressure. We propose the 10K transition occurs in the charge sector, and is a consequence 
of the appearance of the excitonic insulator. Observation of the lattice symmetry breaking 
associated with the excitonic insulator state is an obvious experimental test of our proposal. 
There is a suppression of the spin susceptibility below the transition, and this can arise in 
our theory from the strong coupling between the fermionic charge carriers and the spinons; 
this coupling is given by the large hopping term, t in a t-J model, and we will discuss a 
number of consequences of this term. 

The Z 2 spin liquid state upon which our excitonic insulator state is based has a gap to 
all spin excitations, and this is potentially an issue in applications to k-(ET) 2 Cu2(CN) 3 . 
However, as we have already mentioned above, our theory is easily extend across the critical 
point to the onset of magnetic order. We propose that the spin gap is very small due to 
proximity to this critical point, and that along with the effect of impurities, this can possibly 
explain the large density of states of gapless spin excitations observed in Ref . El Indeed, the T 
dependence of the spin susceptibility fits the predictions of the nearest-neighbor Heisenberg 
antiferromagnet reasonably well, and the latter model is believed to be very close to a 
quantum phase transition at which the Neel order disappears^. 

The central focus of this paper will be on the charged excitations of the Z 2 spin liquid. 
These will spinless fermions of charge ±e: we call the charge +e fermion a 'holon', and the 
charge — e fermion a 'doublon'. We will begin with a theory of the dispersion spectrum of 
single holons and doublons in Section [XTJ The remaining sections of the paper will consider 
various many-body states that can appears in a finite density liquid of these fermions. The 
states we will find are: 

(i) Excitonic condensates: This will be discussed in Section III Here there is an instability 
towards the pairing between a holon and doublon to form a condensate of neutral bosons. 
The only physical symmetry broken by this pairing is the space group symmetry of the 
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triangular lattice, and so the state with an excitonic condensate has a spatial modulation 
in the electronic charge density. We find this state as an instability of the insulating Z 2 
spin liquid, which then reaches an insulating state with an excitonic condensate. However, 
eventually Fermi surfaces of holons and doublons can develop, leading to a metallic state 
with an excitonic condensate. 

(ii) Metals: As the repulsion energy U/t in a Hubbard model becomes smaller, we will find 
that Fermi surfaces of the holons and doublons can appear (as has just been noted in the 
discussion of excitonic condensates). Such states are realizations of 'algebraic charge liquids' 
(ACLs) introduced recently in Ref. [22], and these will be discussed further in Section IV The 
simplest such states inherit the topological order of the Z 2 spin liquid. The areas enclosed by 
the Fermi surfaces in these fractionalized metallic states obey a modified Luttinger theorem, 
whose form will be discussed in Section IV We will also discuss the possible binding of 
holons and doublons to spinons, and the eventual appearance of conventional Fermi liquid 
metallic states. 

(m) Superconductors: The metallic states above are also susceptible to a conventional 
BCS instability to a superconducting state, by the condensation of bosons of charge ±2e. 
This will be discussed in Section |Vj The required attractive interactions between pairs of 
holons (or pairs of doublons) is driven by the exchange of spinon excitations of the Z 2 spin 
liquid. We will compute this interaction, and show that it favors a superconductor in which 
the physical electronic pair operators have a d + id' pairing signature. 



II. CHARGED EXCITATIONS OF THE Z 2 SPIN LIQUID 

As we noted in Section [T| the central actor in our analysis, from which all phases will be 
derived, is the insulating Z 2 spin liquid state on the triangular lattice^^. A review of its 
basic properties, using the perspective of the projective symmetry groupP^I is presented in 
Appendix |Aj 

Here we will extend the S = 1/2 antiferromagnet to the t-J model on the triangular 
lattice and describe the spectrum of the charged holon and doublon excitations of the Z 2 
spin liquid. Because of the spin-charge separation in the parent Z 2 spin liquid, both the 
holons and doublons are spinless, but carry a unit charge of the Z 2 gauge field of the spin 
liquid. The distinct metallic, insulating, and superconducting phases which can appear in a 
system with a finite density of holons and doublons are described in the subsequent sections. 

We consider the usual t-J model on a triangular lattice: 

H = -t^<S a <$ + jY t S i -S j (2.1) 

(ij) ,° (ij) 

where i, j are sites of the triangular lattice. As discussed in Appendix [A] the Z 2 spin liquid 
is realized by expressing the the spins in terms of Schwinger bosons bi a . Here, we introduce 
spinless, charge e, fermionic holons, fi, and spinless, charge — e, fermionic doublons, g^. The 
ha, fu an d g% all carry a unit Z 2 charge. At half-filling, the density of holons must equal the 
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density of doublons, but our formalism also allows consideration of the general case away 
from half-filling with unequal densities. We can now expression the electron annihilation 
operator q q in terms of the these degress of freedom by 

4 = Ml + 9Ua^ (2.2) 

where e a p is the antisymmetric unit tensor. A similar spin-charge decomposition of the 
electron into spinons was discussed in the early study of Zou and Anderson^, but with the 
opposite assignment of statistics i.e. fermionic spinons and bosonic holons. In our case, the 
choice of statistics is dictated by the choice of the Z 2 spin liquid we are considering. 

With the spin liquid background, the charge degrees of freedom acquire a kinetic energy 



from the hopping term. To see this, plug the separation (2.2) into the t-term in equation 



(2.1) 



H t = ~ t E (ftl + ^9}bf) ■ (f}b« + e^ 7 ) + h. c. 

(ij) 

= ~ t E /</X# - 1 E eaf^gUifflrr ~ 1 E ^i^ h \A 

(ij) (ij) (ij) 

-tJ2glfh^ b i b ? + h - c - ( 2 - 3 ) 

(ij) 

Replacing the spinon operators with the mean field expectation values of the spinon pair 
operators discussed in Appendix |Aj we obtain 

/// X:K = - t E fifjBij - t E gtgjBji - t E /;//,. l; ; - t E !i,fj-\r + c. 

(ij) (ij) (ij) (ij) 

= ~ Bt E (f*fi + Ml ) - W E (£9j + 9j9i 

(ij) (ij) 



I E A « (sjf! - iff) - t E A % (fi9j - fj9i 



(ij) (ij) 



where Aij and are the expectation values defined in Eq. (A2). Using Fourier transform 
to convert this to momentum space, we obtain 

H™ F = -2t B E (fkfl + ghk) Re £ fc + 2tA E (ffl/l* + f-k9k) Im & (2.4) 

k k 

where is defined as 

£ k = e ikl + e lk2 + e ik3 (2.5) 

Here we are using the oblique co-ordinate system shown in Fig. [TJ For momentum space, 
we use the dual base with respect to the coordinate system. Therefore k ■ r = k\T\ + k 2 r 2 . 
Momenta (ki,k 2 ), (hi + 2TT,k 2 ) and (ki,k 2 + 2ir) are equivalent. For convenience, a third 
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component of momentum is defined as k% = 



— k\ — &2 to make formulas more symmetric. 




FIG. 1: Oblique coordinates system of triangle lattice^ 



In the complete Hamiltonian, we need a chemical potential \i to couple to the total charge 
density: when this is within the Mott gap, the density of holons and doublons will be equal 
and the system will be half-filled. Also, the double occupancy energy U needs to be included 
since the number of double occupied sites is not conserved. Combining these two terms with 



the kinetic energy (|2.4|), we obtain 
H UF = 



- 2tB (fkft + gbk) Re £ fc + 2tA £ (g\f_ k + f_ k g k ) Im & 

k k 
+ ^Y1 [Ask - fkfk) + U J2 9 k9k 



(2.6) 



The above Hamiltonian can be diagonalized by Bogoliubov transformation. To do so, 
first write the Hamiltonian in the following matrix form: 



(2.7) 



where 



and the kinetic energy matrix 



^k = (gk, f- k ) 



-2tBRe£ k +n+ U 2tAlm£ k 



2tA Im& 



-2tBRe£ k + fi 



The spectrum of quasi-particles is the eigenvalue of the matrix (2.9) 



U U 2 

E k = -2tBReZ k + n+-±d — 



-4t 2 A 2 (Im£ fc ) 



We can redefine the chemical potential \i to absorb the \ term in Eq. (|2.10|). Also, to 



(2.1 



(2.9) 



(2.10) 



make it simpler, define 



He. 



U 



(2.11) 
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Then the dispersion of the two bands looks like 



-2tBRe£ k ± J /I* + (2L4Im£ fc ) 2 



(2.12) 



To study the maximum and minimum position of the about dispersion, we need to first 
identify the stationary points of the dispersion in the Brillouin zone. This can be done 



by noticing that the two k dependent function Re^ and Im^ in equation (2.12) reach 



stationary points simultaneously at the points listed in the first column of Table |T} The 



Steady points 


Re£ fc 


1 Im£fc| 


Et 


± (§*>§*) 
(0,0) 

(7T,7r), (0,vr), (tt,0) 


3 
2 

3 
-1 


§V3 




3tB±y/nl+ (3VstA) 2 
—6tB ± /Xfl 
2tB ± /i fl 



TABLE I: Stationary points of 

maximal and minimal points of the bands must come from the candidates listed in Table |T} 
To determine which one is the maximum or minimum, we compare the energy at those 
stationary points, which is listed in Table [I] along with real and imaginary parts of It is 
easy to see that the position of minimum and maximum depends on the sign of tB and the 
value of yUfi, tA and tB. Let us discuss two extreme case fi^ <C tA,tB and /ifi ^> tA,tB for 
both tB > and tB < 0. 



tB > 

First, consider the extreme case hq 3> tA. In this case the effect of the tA term can 
be ignored comparing to the /xa term in the square root and the shape of dispersion is 
solely determined by the tB lie term, which is listed in the second column of table 
[ij Therefore the maximums locate at ± (|7r, |7r), and the minimal point is (0,0), for 
both the upper and lower bands. 

Then consider the other extreme: /ifl = 0. In this case the energy at the stationary 
points are 



E± = 3tB±3V3\tA\, fc = ± U 7r '^ 7r J 



-6tB, k = (0,0) 
= 2tB, fc=(7r,7r),(0,7r),(7r,0) 



(2.13) 

(2.14) 
(2.15) 



Without the tA term, ( |2.13 ) was the maximum. However, when the 3y/3tA term 
is deducted from it, it becomes smaller and can be smaller than (2.15). Then the 



maximal points of the lower band become (vr,7r), (0,7r) and (vr,0). The requirement 



for (2.15) to be smaller than (2.13) is 



\B\ < V3\A\ 



(2.16) 
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By minimizing the Hamiltonian the mean field parameters A and B we get show that 



the inequality (2.16) holds. Therefore the minimal points are (0, n) and (7r, 0). 



On the other hand, for the upper band the tA term makes equation (2.13) larger so 



that it will never be smaller than (2.14). Consequently the minimal point of upper 



band remains at (0,0). These results are summarized in Table [H 



tB < 

For this case, notice that the dispersions of the bands for parameters tB > and 
tB < are related 

E±(tB) = -ET(-tB) (2.17) 

Therefore for B < the bands are up-side down. Hence the bottom of upper band 
and the top of the lower band are switched (see Table |H|) . 



tB 




Min of e+ 


Max of e 


tB>0 


A*a > \tA\ 


(0,0) 


± 


tB>0 


Atfl < \tA\ 


(0,0) 


(0,0),(0,7T),(7T,0) 


tB < 


Va > Ml 




(0,0) 


tB < 


t*& < \tA\ 


(0,0),(0,7T),(7T,0) 


(0,0) 



TABLE II: Minimal and maximal points of the bands 



According to Eq. (2.17), the dispersion of the bands are flipped when tB changes the 
sign. However, since we are at zero doping, the flip of the bands only switches the name of 
holons and doublons, but the physics is very similar in the two cases. From minimizing the 
spinon mean field Hamiltonian we found that B < 0. Therefore in the following discussion 
we will assume that B < and t > 0. Consequently the maximum of valence band and the 
minimum of electron band comes from the tB < row of Table |TT} when the band gap is 
small /ifj <C \tA\, lower band excitations are around (0,0), and higher band excitations are 



around (0, tt), (tt, 0) and (tt, tt). It is worth noticing that the mixing term in Eq. (2.6), which 
is proportional to tAlm^k, vanishes at these points (see Table |TJ). Therefore, the Bogoliubov 
transformation is trivial at the center of the Fermi sea of the two types of excitations. If we 
assume that the Fermi sea of excitations are small, then the excitations can be identified as 



holons and doublons. According to the definition of quasiparticle (2.8), it carries the same 
electric charge as the original electron. Therefore the quasiparticle excitation in the upper 
band is doublon, and the quasihole excitation in the lower band is holon. Therefore low 
energy effective fields of holon (/) and doublon (h) can be defined as 

/(»o = X) /« e< * r (2 - 18) 

9x(r) = 9k° +g e tq - r , X = 1,2,3 (2.19) 
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where denotes three minimums in doublon dispersion: 

fef = (0, tt), fcf = (vr, 0), fcf = (vr, vr) (2.20) 
The positions of these excitations in the Brillouin zone are indicated in Fig. [2j 



k 2 




FIG. 2: Brillouin zone of the triangular lattice. The first Brillouin zone has the hexagon shape, 
and the two axes are the dual base with respect to the oblique coordinates we are using in Fig. [TJ 
The symbols corresponds to the momenta of minimal energy of excitations. The solid triangles at 
± ^) are spinon excitations (see Appendix |a| ; the solid hexagons at (ir,0), (0, 7r) and (vr,7r) 
are doublons g; the solid square at the origin is the holon /. The hole (bound state of holon 
and spinon) is located at the same position as the spinon (triangles); the electron (bound state of 
doublon and spinon) locates at the solid circles: ± (— §, ), ± (x> — f ) an< ^ ^ (§> §)• ^he bound 
states are discussed in Appendix [Cj 



A. Two-band model 



In this section, we setup a two-band model describing holon and doublon excitations, 
and the Coulomb interaction between them. In the next section, this model will be used to 
investigate exciton condensation in the system when the charge gap becomes small. 

When the density of charge excitations is small, the Fermi pockets of spinon and holon 
are very small in the Brillouin zone. Therefore we can make the approximation that the 
dispersion relation in equation (2.12) takes a quadratic form in the pocket. Without losing 
any generality, we can assume that t > and B < 0. So from Table [TT] holon pocket is at 
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(0,0), and doublon pockets are located at (vr,0), (0, ir) and (vr,7r). Expand the dispersion 
relation in Eq. (2.12) to the quadratic order, we have 



E k ~6t\B\ -/i fl - -t\B\(k 2 x + h 
I 

near the minimum of holon, and 

, , (U 2 A 2 t\B\\ 
E k ~ -2t\B\ +/x fl + ( ^)ki 



(2.21) 



(2.22) 



/i fl 2 

near the minimum (0, 7r) of doublon. The dispersion relation near the other two doublon 



minimums can be obtained by simply rotating Eq. (2.22) by 120° degrees. 



From Eq. (2.21 ) and (2.22) the gap between the two bands can be calculated by subtract- 



ing the maximum of holon band from the minimum of doublon band 



^ind = 2/ifl — 8t\B\ 



(2.23) 



and it is an indirect gap because the maximum of holon band is at (0, 0), while the minimums 
of doublon bands are at other places: (0,7r), (tt,0) and (7r,7r). If we do not include any 
interactions, the system consists two bands: holon band and doublon band. If the gap 
£/ in d is positive, it is a semiconductor with an indirect gap. If the gap is negative, it is a 
semimetal. 



From Eq. (2.21 ), we see that the dispersion relation of holon is isotropic with the effective 
mass: 

"* = Mb\ (2 - 24) 



From Eq. (2.22) we see that the dispersion of doublon is not isotropic. For the pocket at 



(7r, 0), the effective masses in x and y directions are 



m dx 



16t 2 A 2 



t\B\ 



/'ft 



m dy 



3t\B\ 



(2.25) 
(2.26) 



Note, that m dy = m^. When U in d is tuned around zero, the Coulomb interaction will lead to 
exciton instability. The change of ?7; n d is small comparing to t, so we can assume the masses 
are constant in our problem. At the point [/i n( j = 0, we have n$ — At\B\. Therefore 



m dx 



m dy 
m dx 



1 

3 



(2.27) 



In the spin liquid phase we have A > \B\, therefore the above ratio is larger than 1. 

The dispersion relation for the other two doublon pockets can be obtaned by rotating the 



10 



pocket by ±120°. 

The exciton problem has been well studied in the context of semimetal to semiconductor 
transition^, where people usually use electron operators for both the valence and conducting 
band. Therefore here we want to choose our notations to be consistent with them. To achieve 
that, we just need to switch the creation and annihilation operator for holon. To be clear, 
we use h and d for new holon and doublon operators respectively. 

h k = ft, 4 = 9k (2.28) 
For the doublon bands, we measure the momentum from the minimum in the pockets: 

dlk = dk+( n ,0), dih = djfc+(0,7r)> dsk = dk+(n,ir) (2.29) 

With the help of those effective masses, the dispersion relation for holon and doublon can 
be expressed. 

k 2 

e h {k) = U ind - — (2.30) 
2m h 

k 2 k 2 

e d2 (k) = 7r ^ + 7 ^ L - (2.31) 
2m dx 2m dy 

/IN (^ cos 120° - A;„ sin 120°) 2 (k x sin 120° + k v cos 120°) 2 

e d i(k) = K — y - '- + y — y - '- 2.32 

2m dx 2m dy 

fu , _ {k x cos 120° + k y sin 120°) 2 {-k x sin 120° + k y cos 120°) 2 

£d3{k) — 1 (2.33) 

2m dx 2m dy 

Therefore the kinetic energy of the system is 

#kin = e h{k)h\h k + e dX {k)d{ k d Xk (2.34) 

The exciton instability is because of the Coulomb interaction between doublon and holon. 
In general, the Coulomb interaction is 

1 e 2 

Hmt = 2 \ r . _ r A* C ia C )p C iP (2.35) 

where the electron density can be expressed by holon and doublon operators as the following 

J2 c l c ^ = h % + 4di (2.36) 
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Therefore Eq. (2.35) becomes 

o2 



#int = \Y. \ r e _ r .\ + d ldi) ■ (h)h 3 + d)d^j (2.37) 

= Yl ~T ( h l h k+Q + d Uk+q) ■ (hk'hk'-q + 4'4'- g ) (2.38) 
kk'q Q 

To analyze this interaction, we can use the screened Hartree-Fock approximation^ For 
the Hartree term, the unscreened Coulomb interaction is used; for the exchange term, the 
screening from other electrons is included. This approximation results in the following 
functional of interaction energy 



H ^ = E ^ + (4dk +q )) ■ + (4'4'- ? )) (2.39) 

Aire 2 
e(q)q 



4 ^ 

H ex = - (( h k h k'-q) ( h i h k+ q ) + (dl d k'- q ) (dl'dk+q 



k,k' ,q,aP 

+ (hldy- q ) (dlhk+q) + (dlhk^q) (hld k+q )) (2.40) 



where (■) denotes the average under the mean field state. 



The exciton is formed by pairing of holon and doublon. In the interactions in Eq. (2.39) 



and (2.40), only the last two terms in exchange interaction are relevant for holon and doublon 
pairing. Therefore in the exciton problem we can include only these terms and neglect all 
the others. This is called dominant term approximation^. 

4 ^ 

#dta = -J2 7^2 (( h ldk>- q ) (4,h k+q ) + {d\h k ,- q ) (hld k+q )) (2.41) 

kk'q 



Since the interaction is invariant in parity transformation, equation (2.42) can be further 
simplified as 

4 ^ 

H DTA = -2 7(^-2 (hUk'-q) (dl,h k+q ) (2.42) 

kk'q WW 

The interaction can be further expanded near the minimums of holon and doublon pockets 
as the following 

HdTA = ~ E e(q + 4)(q + fcf) 2 ( h * d ™-«) (2-43) 

kk'q,\ KH ^ A ' 

In this equation k, k' and q are both small momenta. The q dependence in the coefficient 
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of this interaction can therefore be neglected: 



8 - 2 (2.44) 



and we define this coefficient as g. The interaction in equation (2.43) can then be expressed 
as 

kk'q,\ 

III. EXCITON CONDENSATION 

In the previous section we discussed the dispersion relations of holon and doublon ex- 
citations and the interaction between them. If momentum is measured from the center of 
the pockets, the holon pocket is a sphere centered at (0,0), while the doublon pockets are 
three ellipses centered at (0, 0), long axes of which are 120° to each other. Since the system 
is neutral, the area of the holon pocket is three times larger than the area of each doublon 
pocket. However, the long axis of doublon pockets can still be longer than the radius of 
holon pocket if the aspect ratio of doublon pockets is greater than 3:1. Therefore there can 
still be pairing at the places where the two Fermi surfaces are close. Between holon and one 
type of doublon the pairing is only good in one direction, but if the holon is paired with 
both three types of doublon, there is a good pairing in three directions, which covers a large 
portion of the Fermi surface. We will see that this mismatch in Fermi surfaces lead to a 
richer phase diagram compared to the normal-superconducting problem. 

Since all the pockets are centered at (0, 0) (when momentum is measured from the min- 
imums of pockets), we assume that the exciton pairs also have total momentum (0,0). In 
other words, it is not a LLFO state in which smaller pocket is shifted toward the edge of 
larger pocket. Therefore we have three pairing order parameters between holon and each 
type of doublon 



Note, that the pairing is s-wave because the interaction (2.45) only has s-wave component. 
According to this order parameter, the mean field decomposition of the interaction (2.45) 

is 1 

h dta = Ys A *Yl 4a + h - c - + - E i a ^i 2 ( 3 - 2 ) 

\ 7„ 9 \ 



Combining this with the kinetic energy in Eq. (2.34) and the chemical potential, the full 
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Hamiltonian is 



H =-Hkin + -f^DTA — A 1 



Afc — iV( 

^[e fc (fe) - /tx]4^fc + ^2[td\{k) - ^}d\ k d Xk 

k k,X 

+ J2 Ax J2 + h. C. + -J2 I a aI 2 + VN 



9 



(3.3) 



A k X 

This Hamiltonian can be diagonalized using the Bogoliubov transformation. First, write 



Eq. (3.3) in the matrix form: 

H = *lM k *k + - ^ | A A | 2 + fiN 



where 



^fc — (hk dik d 2 k d 3 k) T 



A3 \ 



(e h (k)-n A* A* 

Ax e dl (k)-fx 

A 2 e d2 {k)-n 

V A 3 e d3 (k)-fij 

A Bogoliubov transformation is then a unitary transformation 



(3.4) 



(3.5) 



(3.6) 



(3.7) 



that diagonalizes the matrix Mk- With the new quasiparticle operators the Hamiltonian 
is diagonalized as 



H = J2Y1 E ik(*lk a *k + " E |A A | 2 + fj,N 

j„ — 1 9 \ 



(3. 



k i=l 



where are the four components of ^k, and Eik are the four eigenvalues of the matrix Mk- 
Therefore in the ground state, the state of quasiparticle is filled if < 0, or empty if 
E ik > 0. The ground state energy is 



(H)= £ ^ + -^|A A | 2 + /iiV 

rr ^-n 9 \ 



(3.9) 



E ik <0 



Becuase the three doublon pockets are identical, we assume that the absolute values of 
the order parameters are the same for three pockets. Therefore the pairing order parameters 
can be characterized by one absolute value and three phases. 



A> = Ae 



(3.10) 
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After we have made the approximation that only excitations near the minimums of the 



pockets are considered, the kinetic energy (2.34) and interaction energy (2.45) has the {7(1) x 



{7(1) x {7(1) symmetry for changing the phases of three doublons independently. Therefore 
the mean field energy (3.3) does not depend on the phases of order parameters 6\. In other 



words, for this simplified Hamiltonian, a state with order parameter in Eq. (3.1) breaks 



the {7(1) x {7(1) x {7(1) symmetry spontaneously, and the states with different phases are 
degenerate. We will see in the next section that this degeneracy is lifted once we add the 
other parts of Coulomb interaction back. 



With this {7(1) x U(l) x {7(1) degeneracy, the ground state energy in Eq. (3.9) only 



depends on the absolute value of A. The ground state is obtained by minimizing the ground 
state energy with respect to order parameter A. The results are described below 

When [/md is large and positive, the system is in an insulating phase (semiconductor 
phase), with no Fermi surface and A = 0. Then as the gap decreases, the system goes 
through a second-order phase transition at a critical value U® nd and enters a insulating phase 
with no Fermi surface but A ^ 0. In this phase the binding energy of exciton overcomes 
the charge gap, so there are doublon and holon excitations. However, the excitations near 
the Fermi surface bind into exciton, and the gap is large comparing to the Fermi energy 
such that the Fermi surface is fully gapped although there is a mismatch between holon and 
doublon Fermi surfaces. This state is an insulating phase because the charged quasiparticle 
is gapped. We call this phase an excitonic insulator. Then if the gap is further decreased, 
the system goes through a second-order phase transition at some value U^ d < into a phase 
still with A / 0, but there are Fermi surfaces in the spectrum of quasiparticle excitation. 
Hence this is a conducting phase with A ^ 0. A typical image of the Fermi surfaces in the 
system is shown in Fig. [5} Eventually, the system goes through another first order phase 
transition at U^ nd < U? ad into a state with A = 0. This is a simple semimetal conducting 
state. A typical phase diagram for given masses and g is shown in Fig. [3j For given masses, 
the two dimensional phase diagram for different g and U- m d is shown in Fig. |4} 

In phase with the excitonic condensate, the ground state has a {7(1) x {7(1) x {7(1) 
degeneracy because of the corresponding symmetry in the simplified Hamiltonian. This 
degeneracy will be lifted to a {7(1) degeneracy once we take into account the parts in the 



interaction (2.38) which we ignored in the dominant term approximation. In other words, 



the full interaction breaks {7(1) x U (1) symmetry, and only one {7(1) is spontaneously broken. 



Besides the dominant term in Eq. (2.41), the next leading term is the direct interaction 



in Eq. (|2.39|), because it is unscreened. That term can be expressed using charge density 

(3.11) 



operator as 



H. 



dir 



J2 (pi) ■ (p-i) 



q^0 



where p q is the density operator 



Pg 



(3.12) 
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FIG. 3: Phase diagram at fixed masses and interaction g. This plot shows the exciton pairing gap 
A and the area of hole-like quasiparticle Fermi surface as a function of the charge gap Ui n d- It is 
plotted at nth = m^y = Sm^x, and g/rrih = 0.007. J7j n( j and A are measured in unit of the energy 
scale in the system: where A is the cutoff momentum; the Fermi surface area is in unit of it A 2 . 
Phase I, II, III and IV in the diagram are normal insulating phase, excitonic insulator, excitonic 
metal and normal conducting phase respectively. 



In a state with translational invariance, the charge density is a constant, and p q vanishes 
at any q 7^ 0. However, with an excitonic condensate, translational symmetry is actually 
broken. This can be understood by the following ways: First, in the excitonic condensate, 
there is pairing between holon and three types of doublon. The total momenta of these pairs 
are not zero: a pair of one holon and one doublon d\ has total momentum k^. Therefore 
the translational symmetry of the system is broken. Also, one can look at the change of 
Brillouin zone. Originally the doublon pockets and holon pocket are at different places in the 



Brillouin zone. However, after the Bogoliubov transformation (3.7), the holon and doublon 



modes are mixed together. The only way this is allowed is that the Brillouin zone shrinks 
to one quarter of its original size, so that are all reciprocal vectors. The shrinking of 
Brillouin zone means that the unit cell in real space is enlarged, and the original translational 
symmetry is broken. 

Because of the broken symmetry, we can look for non-vanishing expectation value of 
density operator at q = k^. Let us consider fcf = (tt, 0) as an example. Again we consider 
only the excitation near the center of pockets. The expectation value of density operator is 

(p<t,o)) = E ((4<*3*) + (4kd2k)) (3.13) 
k 

In our ansatz we do not have this average put in explicitly. However, because the order 
parameter A\ already breaks the U(l) x (7(1) X U(l), we can get a non-vanishing value for 
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FIG. 4: Phase diagram at fixed masses. It is plotted at nih = may = Srridx- Phase I, II, III and 
IV in the diagram are normal insulating phase, excitonic insulator, excitonic metal, and normal 
conducting phase respectively. 




FIG. 5: Fermi surface in the excitonic metal. Region I is the hole-like Fermi surface of the quasi- 
particle, and region II is the particle-like Fermi surface of the quasiparticle. Because the system 
is charge neutral, the area of region I is equal to the area of region II. Note, that this diagram is 
drawn in the reduced Brillouin zone. In presence of the pairing, the holon and doublon modes are 
mixed, and the Brillouin zone shrinks to 1/4 of its original size such that the holon and doublon 
pockets locate at the same momenta. 
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d\ k dsk ) from the ansatz. 



In the mean field state, the average in Eq. (3.13) can be calculated using the Green's 



function of doublon. Consider the Green's function matrix 



Then according to the Hamiltonian in Eq. (3.4), we have 



(3.14) 



(3.15) 



Then the average in Eq. (3.13) can be obtained by taking the inverse of the above matrix 

uj -E dl (k) + n 



d 2k d 3k 



du 

2vr (w - E kl )(u - E k2 )(uj - E k3 )(u - E k4 



-A*A 3 oc A^Ag 



(3.16) 



Therefore the average of density operator in equation (3.13) is 



<p K0 )> oc Re A^Aa oc A 2 cos(0 2 - 3 ) 



(3.17) 



The expectation value of charge density operator at the other two momenta can be evaluated 
similarly. Finally we get 



4e_ 2 

7T 



H dir oc — A 2 [cos 2 (0! - 2 ) + cos 2 (6 2 - 6 3 ) + cos 2 (0 3 



Oi)] 



(3.18) 



It can be proved that this energy has a minimum at 

2n 



5tt 



03 - 02 = 02 - 0! = — , Or - 



(3.19) 



Such a state breaks the time reversal and parity symmetry because of the complex phase, 
but the rotational symmetry of the original lattice is preserved. Finally from the density 
operator one can conclude that there is a charge order in the real space shown in Fig. |6j 
From the figure one can see that the unit cell is a 2 x 2 plaque. This is consistent with the 
conclusion that the Brillouin zone shrinks to one quarter of its original size. 



IV. FRACTION ALIZED METALS WITH Z 2 TOPOLOGICAL ORDER AND LUT- 
TINGER RELATIONS 

The simplest route from the insulating Z 2 spin liquid to a metallic state is by formation 
of Fermi surfaces of the holons and doublons. We can easily detect such a transition as a 
function of charge density and system parameters by examining the dispersion relations of 



the charged excitations in Eq. (2.12). These Fermi surfaces have quasiparticles of charge ±e, 



but which are spinless. So the resulting state is not a conventional Fermi liquid, but a Z 2 
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FIG. 6: Charge order pattern, t/n can be either positive or negative depends on if the relative phase 
in equation (3.19) is 2 " 



or 



'algebraic charge liquid^. The Z2 topological order of the insulator survives into the metal, 
and quasiparticles carry Z 2 gauge charge. This section will address two related questions: 
what are the Luttinger relations which constrain the areas enclosed by the Fermi surfaces 
of this fractionalized metal, and how does this metal evolve into a Fermi liquid ? 

We will examine a route to the Fermi liquid here which is similar to that considered 
recently by Kaul et al^ for related states with U(l) topological order on the square lat- 
tice. We allow for the binding between the holon and doublon quasiparticles described in 



Section |II A[ and the bosonic spinons described in Appendix |A 1| This binding is mediated 
not by the gauge force (which is weak in the fractionalized phases we are discussing here), 



but by the hopping matrix element t in Eq. (2.3). This matrix element is larger than J, 



and so the binding force is strong. A computation of the binding force for the low energy 
holons/doublons and the spinons is presented in Appendix [Cj and position of these bound 
states in the Brillouin zone is illustrated in Fig. [2] The resulting bound states will have the 



quantum numbers of the electron operator in Eq. (2.2). Fermi surfaces of these 'molecular 



bound states can also formal, with electron-like quasiparticles which carry charge ±e and 
spin 1/2. Indeed it is possible that these bound states have a lower energy than the holons 
or doublons, in which case the first metallic state state will not be an algebraic charge liquid 
(with holon/doublon Fermi surfaces), but have electron/hole Fermi surfaces leading to a 
'fractionalized Fermi liquid^^SEH ( see below and Fig. [7]). 

So as in Ref. [221 the fractionalized metallic states being considered here can have a variety 
of distinct Fermi surfaces. As in a usual Fermi liquid, each Fermi surface can be considered 
to be a surface of either charge e or charge — e quasiparticles, with the two perspectives 
related by a particle-hole transformation. There is also a corresponding change in the region 
of the Brillouin zone which is considered to be 'inside' the Fermi surface, and the value of 
the area enclosed by the Fermi surface. For convenience, let us use a perspective in which 
all Fermi surfaces are described in terms of charge — e quasiparticles, the same charge as an 



electron. Thus in the notation of Eqs. (2.28) and (2.29), we are considering Fermi surfaces of 



the h, d, and c a fermions. Let us by Ah, Ad, and A c the area of the Brillouin zone enclosed 
by these Fermi surfaces, divided by a phase space factor of 47r 2 . We are interested here in 
the Luttinger relations which are satisfied by these areas. 
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As has been discussed in Refs. 25 28 29], there is a one Luttinger relation associated with 
each conservation law, local or global. We have the conservation of total electronic charge 
which we can write as 

y ^ (.</, </; fif) •<• (4.1) 

s i 

where N s is the number of sites, and x is the density of doped electrons away from half-filling 
i.e. the total density of electrons in the band is 1 + x. There is also the local constraint 



associated with the Schwinger boson representation in Eq. (2.2) 



f}fi + gl9i + J2 b l b ? 



(4.2) 



A key observation is that in any state with Z 2 topological order, the global conservation 



law associated with Eq. (4.2) is 'broken', due to the appearance of the anomalous pairing 



condensate (eapbfbj). Consequently^, the Luttinger relation derived from Eq. (4.2) does 
not apply in any such state. 

We can now derive the constraints on the Fermi surface areas placed by the above relations 
using an analysis with exactly parallels that presented in Ref. |25j We will not present any 
details here, but state that in the final result we simply count the contribution of each term in 
the conservation law by counting the number of all such particles inside every Fermi surface, 
whether as bare particles or as components of a molecular state. It is important in this 
counting to also include the number of bosons i.e. the Luttinger relation applies to bosons 
too™M of course for the boson number counting, there is a contribution only from the 
fermionic molecular states. From such an analysis, after properly performing the particle- 
hole transformations associated with the mapping to the charge — e h and d particles, we 



obtain from Eq. (4.1) the relation 



2 A c + A + A d = 1 + x 



(4.3) 



This relation applies in all non-superconducting phases. The corresponding relation from 



Eq. (4.2) is 



Ah + Ad = 



(4.4) 



We reiterate that Eq. (4.4) applies only in a state without Z 2 topological order. Indeed, a 
state without Z 2 order cannot have Fermi surfaces of fractionalized excitations, and so the 
only sensible solutions of Eq. (4.4), when it applies, are Ah — A g — and Ah — A„ — 1. 



With these Luttinger relations at hand, we can now describe their implications for the 
various phases. These phases are illustrated in Fig. [7j and described in turn below. 

a. Z 2 insulator: This is the state reviewed in Appendix |A} The chemical potential is 
within the band gap of the dispersion in Eq. (2.12), as shown in Fig. [7J We have x = 0, 
Ah = 1, Ad = 0, and A c = 0, and so the relation (4.3) is obeyed. This relation also applies 



when an excitonic condensate appears: the only change is that there is mixing between the 
h and d bands. 
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Z 2 insulator 



doublon metal 




doublon-electron metal Fractionalized 

Fermi liquid 




holon- electron metal Fractionalized 

Fermi liquid 



Fermi liquid 

FIG. 7: Schematic illustration of the spectrum of fermionic holon (h), doublon (d), and electron 
(c Q ) excitations in various phases which descend from the Z2 spin liquid. The band structures are 
merely represented by cartoon sketches with the purposes only of indicating the filled states and 
the location of the chemical potential. More realistically, the Fermi surfaces will form around the 
respective dispersion minima in Fig. [2] 
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b. Z 2 metals: A variety of metallic phases can appear, which are all algebraic charge 
liquids^, depending upon the presence/absence of h, d, and e Fermi surfaces, as shown in 
Fig. [7} All these phases have Z 2 topological order, and so only the sum-rule (4.3) constrains 



the areas of the Fermi surfaces. As an example, we can dope the Z 2 insulator with negative 
charge carriers (x > 0). This will lead to a 'doublon metal' with Ah = 1, Ad = x and A c = 0. 
As the energy of the bound state between a doublon and spinon is lowered (with increasing 
t), a Fermi surface of electrons can appear, leading to a state with Ah — 1, < Ad < x 
and < A c < x/2. This fractionalized state with both h and e Fermi surfaces (a 'doublon- 
electron' metal) is the particle-hole conjugate of the 'holon-hole' metal discussed recently in 
Ref. [221 The location of the electron Fermi surface in this metal is computed in Appendix [C| 
and will reside at pockets around the electron dispersion minimum in Fig. [2j 

c. Fractionalized Fermi liquid: With a continuing increase of the doublon-spinon bind- 
ing energy in the doublon-electron metal just discussed, there is a transfer of states from 
within the doublon Fermi surface to within the electron Fermi surface, until we reach the 



limiting case where Ad = 0. At the same time we have Ah = 1, and then Eq. (4.4) requires 
A e = x/2. So the chemical potential is in the gap between the doublon and ho Ion bands, as 
shown in Fig. [7j This state is the fractionalized Fermi liquid of Refs. l2"o r f2"Tl 

d. Transition to the Fermi liquid We now continue decreasing the energy of the electron 
bound state, in the hopes of eventually reaching a conventional Fermi liquid. One plausible 
route is that illustrated in Fig. [7j As the electron energy is lowered from the fractionalized 
Fermi liquid, the chemical potential goes within the lower holon band, and we obtain a 
fractionalized holon-electron metal. Note that this state has both an electron-like Fermi 
surface with charge — e carriers (as measured by the contribution to the Hall conductivity), 
and a holon Fermi surface with charge +e carriers. Eventually, the chemical potential goes 
below both the holon and doublon bands (see Fig[7J), and we obtain another fractionalized 
Fermi liquid with Ah = Ad = and A c = (1 + x)/2. Note that this is the same area of an 
electron-like Fermi surface as in a conventional Fermi liquid. However, the Z2 topological 
order is still present, which is why we denote it 'fractionalized': the fractionalized holon and 
spinon excitations are present, but all of them have an energy gap. Eventually, however, 
we expect a confinement transition in the Z 2 gauge theory, leading to a conventional Fermi 



liquid with no fractionalized excitatons. Because the Luttinger relation in (4.4) is saturated 
by the electron contribution, we expect this confinement transition will be in the universality 
class of the 'even' Z 2 gauge theory, rather than the 'odd' Z 2 gauge theory^^l 



V. SUPERCONDUCTING STATES 

The projected symmetry group analysis can be extended to analyze the pairing symmetry 
in possible superconducting states. In a superconducting state, the pairing symmetry can 
be determined by the relative phases of the expectation value 



Ay = (eapCiaCjp) (5.1) 
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at different directions. On the other hand, the PSG determines how the operator Cj CT trans- 



forms in rotation, and consequently determines how Eq. (5.1) transforms when the bond (ij) 
is rotated to another direction. Therefore the PSG of the electron operators can determine 
the pairing symmetry. 

In our Schwinger-boson-slave-fermion ansatz, the fermionic holon and doublon excitations 
form Fermi seas around the bottom of their bands. If there are appropriate attractive forces, 
holons and doublons will become superconducting, and the pairing of holons or doublons 



will make Eq. (5.1) nonzero because the spinon is already paired. Similarly, the pairing 



symmetry of holon and doublon can also be determined by their PSGs, which are related 



to the PSG of electron according to Eq. (2.2). Therefore the pairing symmetry of electrons 
and quasiparticles are related. 

When determining the PSG of holon and doublon operators, it is assumed that the 
physical electron operator Cj<j is invariant under symmetry transformations. However, in the 
superconducting state the gauge symmetry of electron is broken and the electron operator 
may need a gauge transform after a lattice symmetry operation to recover. The pairing 
pattern is directly related to the gauge transform associated with rotation. Assume that the 
electron operator transforms in the following way after a Ril operation: 



c i(T e ia (5.2) 



Then after the rotation the superconducting order parameter (5.1) becomes 



Ajj -> A ije 2ia (5.3) 

This equation shows the pairing type of the superconducting state. 

If we want a superconducting state with translational and rotational symmetries, the 



order parameter (5.1) should recover after acting Rtl three times. This requires the phase 
a satisfies 

(e 2ia ) 3 = 1 (5.4) 

which implies that a = 0, |, It is easy to see that these choices correspond to s-wave 
and d x 2_ y 2 ± id xy pairing patterns. (Note, that the <i-wave pairing breaks time reversal 
symmetry) . 



On the other hand, Eq. (5.2) also determines the PSG of holons and doublons, and 



therefore determines their pairing patterns. For the d — id pairing, the result is in Table III 



For the s-wave pairing, the result is in Table IV 



From Table III and IV we can see that p-wave pairing of holon and doublon will result 
in d-wave pairing of electron; /-wave pairing of holon and doublon will result in s-wave 
pairing of electron. In the rest of this section we will show that the interaction mediated by 
spinon excitation leads to a p-wave pairing symmetry for holon and doublon, and results in 
a (i-wave pairing for electron. 

In the spin liquid state, the excitation of spin degree of freedom is described as Schwinger 
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TABLE III: PSG and pairing symmetry of charge excitations. 
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TABLE IV: PSG and pairing symmetry of charge excitations. 



boson. The holon acquires a kinetic energy from the original hopping term by replacing 
the product of spinon operators by the mean field average. Moreover, beyond the mean 
field theory this term also describes interaction between holon and spinon excitations. This 
type of interaction results in an interaction between two holons by exchanging a pair of 
spinons. Detailed calculation will show that this interaction is attractive and can lead to 
superconducting pairing of fermionic holons. 

In Appendix [Bj we calculate the scattering amplitude between two holons by exchanging 
two spinons shown in Fig. [8j Also, a contribution from 0(1/N) fluctuation of spinon mean 
field solution is included to make the result consistent with the constraint that there is one 
spinon on every site. 




FIG. 8: Interaction 



The scattering between two holons by exchanging two spinons can be represented as a 
Feynman diagram in Fig. [8] The vertex in this figure comes from the hopping term in 
equation (2.3). For Cooper pairing instability, we are interested in the following situation: 
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the incoming and outgoing momenta are on the Fermi surface; the total incoming and 
outgoing momenta are zero; the energies of incoming and outgoing particles are on shell. 
Therefore the interaction T is a function of the angle 9 between incoming and outgoing 
fermions. In Appendix |Bj we calculate the scattering amplitude of this Feynman diagram. 
Also, a contribution from 0(1 /N) fluctuation of spinon mean field solution is included to 
make the result consistent with the constraint that there is one spinon on every site. The 
result is the following 



When m s ^> q, we get 



When m s <C q, we get 



t 2 ( „<f 



r~- T Am s + C^) (5.5) 



t 2 



^ (5-6) 



where A is a UV cutoff of the system, m s is the spinon gap, and q = 2fcpsin| is the 
momentum transfered in the interaction. 

Since the incoming and outgoing particles are fermions, we need to antisymmetrized T: 

f(0) = r(0) -T(-0) (5.7) 



The antisymmetrized interaction in the limit of equation (5.5) is proportional to cos0, and 



therefore has only p-wave component. The antisymmetrized interaction for the other limit 



(Eq. (5.6)) is a mixture of odd-parity wave components, but the p-wave component is dom- 
inant. Consequently this interaction will lead to a p-wave superconducting instability for 
holon, and result in a <i-wave pairing superconducting state in terms of electron according 
to table ED 



VI. CONCLUSIONS 

This paper has presented a route for the transition from an insulator to a metal for 
electrons on a triangular lattice. Our primary purpose was to pay close attention to the an- 
tiferromagnetic exchange interactions, and the consequent resonating valence bond (RVB) 
correlations in the insulator. We were particular interested in the fate of these RVB cor- 
relations in the metal. Our approach borrows from a recent theory^ of 'algebraic charge 
liquids' (ACL), and can be straightforwardly generalized to other lattices in both two and 
three dimensions. As we noted in Section [TJ our theory of the insulator-to-metal transition is 
complementary to the DMFT theory^ of the metal-to-insulator transition. The latter begins 
with a theory of the correlated Fermi liquid, but pays scant attention to RVB correlations, 
which we believe are particularly important on the triangular lattice. 

We have introduced a variety of intermediate phases in the route from the insulator to 
the metal, and so we now present a brief review of their basic characteristics. 

First, let us discuss the possible phases of the insulator and their charged excitations. 



25 



We used a Z 2 spin liquid stattP 131 with bosonic spinons as a central point of reference. 
This state was reviewed in Appendix [Aj and its low energy excitations are gapped S —1/2 
neutral spinons z a which are located (in our gauge choice) at points in the Brillouin zone 
shown in Fig. [2} Its charged excitations consist of spinless, charge +e ho Ions /, and charge 
— e doublons g, which were described in Section [XTJ their positions in the Brillouin zone 
(in our gauge choice) are also in Fig. [2j We also discussed, Appendix [C] the binding of 
holons/doublons with spinons to form charged, 5* = 1/2 excitations which are hole/electron- 
like; the location of these bound states in the Brillouin zone is indicated in Fig [2} This 
binding binding between gapped excitations is caused by the t term in the t-J model, and 
so the attractive force is quite strong. Consequently, it may well be that the lowest charged 
excitations of a Z 2 spin liquid are electrons and holes, and not doublons and holons. Note 
that the binding is present already in the deconfmed spin-liquid state, and is quite distinct 
from the confinement transition. 

Although, we did not say much about this explicitly, our analysis of excitations in the 
Z 2 spin liquid is easily continued into confining insulators with Neel or VBS order. The 
transition into the Neel state is driven by the condensation of the z a spinons, as described in 
Ref. [H The charged excitations in the Neel state are simple continuations of the holons and 
doublons in the spin liquid, as has been discussed in some detail previously for the square 
lattice^. The charged excitations in the VBS state will be the electrons and holes described 
in Appendix [Cj 

Continuing our discussion of insulating states, our primary new result was the prediction 
of an insulator with an excitonic condensate in Section IIHI We described the excitonic 
instability starting from the Z 2 spin liquid, mainly because this led to the simplest technical 
analysis. However, we expect very similar excitonic instabilities to appear also from the Neel 
and VBS insulators. One physical manifestation of the excitonic instability was a breaking 
of the space group symmetry of triangular lattice. In our mean-field analysis we found the 
translational symmetry breaking depicted in Fig. [6j which serves as a test of our theory. 

Moving onto metallic states, we found a plethora of possibilities depending upon precisely 
which instabilities of the insulator was present (Neel, VBS, excitonic), and which of the 
charged excitations of the insulator had the lowest energy (holon, doublon, hole, electron). 
Some of these states were algebraic charge liquids, while others were Fermi liquids. However, 
the Fermi surface areas in all these states satisfied either the conventional Luttinger relation, 



or a modified Luttinger relations, and these were discussed in Section IV 



Finally, we also considered the onset of superconductivity due to pairing of like-charge 
fermions in Section [V] Again, for technical reasons, it was simplest to describe this instability 
using the holon/doublon excitations of the Z 2 spin liquid, but similar results are expected 
from instabilities of the other states. 

The possible experimental applications of our results to organic compounds was already 
discussed in Section [j] and so we will not repeat the discussion here. Here we only reiterate 
that our model for k-(ET) 2 Cu2(CN) 3 is a Z 2 spin liquid, close or across its transition to a 
state with Neel order, with an excitonic instability at low temperatures. 
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APPENDIX A: Z 2 SPIN LIQUID 

Consider the model on a triangle lattice with anti-ferromagnetic interaction between 
nearest neighbor sites. 

H = JJ2 S *- S i ( A1 ) 

m 

We use the coordinates system labeled in Fig. [TJ and the Schwinger-boson mean field ansatz 
with the following ground state expectation values: 



b l b i° ) ' A * = ( € <*Mp \ (A2) 



where bi a is Schwinger boson operator. 



2 

a/3 



The ansatz (A2) stays the same if both the spinon operators and A and B parameters 



are transformed under the following gauge transformation 



u 7.rr ? & Vr. 



id 



Aij -> e'^-^Aij (A4) 

Because of this gauge freedom, ansatz with different A^ and B^ configurations may be 
equivalent after a certain gauge transform. Therefore the symmetry of A^ and B^ fields 
may be smaller than the symmetry of the ansatz. In particular, we consider the symmetry 
group of the triangle lattice. (We only consider ansatz that does not break the lattice 
symmetry). When a group element g acts on the ansatz, an associated gauge transform 
G g may be required for the A^ and B^ fields to recover. This gauge transform forms the 



PSG. According to equation (A4), gauge transform G can be defined using the phases 4>i{G). 
Consequently the PSG can also be specified using phases of the gauge transformations. The 
phases of the gauge transform associated with symmetry group element g is denoted by (fi g . 

For triangle lattice, the lattice symmetry group can be generated by four elements: T\ 
and T2, which are the translations along the 1 and 2 direction; Rjl, which is the rotation of 
|; a, which is the reflection operation respect to the r 1 = r 2 axis. In the coordinates system 
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used in this note, these elements can be represented as 



Ti : (ri,r 2 ) -> (ri + l,r 2 ) 
T 2 : (ri,r 2 ) -> (r 1; r 2 + 1) 
cr : (ri,r 2 ) -> (r 2 ,ri) 
-Rf : (j"i,r 2 ) -> (ri - r 2 ,r 1 ) 

It can be proved that there are two possibilities for PSG of 2D triangle lattice. We are 
interested in one of them, which is called "zero-flux" ansatz. The zero-flux PSG can be 
described as 

71 

<hi = 0T2 = <\>a = 0, (f)R n/3 = ~ (A6) 

These are also shown in Table IVl 

Assume that the system is slightly doped with holes with hole concentration 5^0. 
Hence the spinon state is not changed. We describe the system by introducing fermionic 
holon 

Ci„ = h\b ia (A7) 
where Cj CT is the electron and hi is holon operator. To preserve the physical electron operator, 



the holon needs to transform as follows in the gauge transformation (A4) 

hi -> e irt>i hi (A8) 
The PSG of holon is summarized in the second row of table |V] 



Operator/Field 


T x 


T 2 


a 


Re. 

9. 


ha 
fi 


fi 


ha 
fi 


ha 

fi 


iha 
ifi 



TABLE V: PSG 



The PSG of spinon operators (A6) determines how and transforms under lattice 
symmetry operations, and therefore reduces the degree of freedom of parameters in the 
ansatz to only two free ones. <ft tl = <Pt 2 = implies that the ansatz is translational invariant. 
By definition and Bij has the following properties 

Aij = -A j{ , B^ = B) t (A9) 

Therefore <p a = and <pR n/3 determines the structure of the ansatz: B^ = B is the same 
real number on every bond; A^ = ±A, and the sign is shown in Fig. [9] 

1. Low energy theory of spinons 

In this section we analyze low energy excitations of spinon. 
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FIG. 9: Configuration of A^ in zero- flux stateP. The arrows on the bonds show the direction 
along which > 0. 



For spinon Hamiltonian, if we take the mean field decomposition of Eq. (Al) using pa- 



rameters (A2), the result is 



Hmf =JJ2 ( A >, A >i + B *iAj + h - c) 

(ij) 

+ J E (Aj a u - B h B a + h - c) - a* E ( E 6 U» " K ) 



(A10) 



where 



A* = E ^ = E e «A« 6 . 



(All) 



a/3 



Plug in the "zero-flux" type of ansatz and write the Hamiltonian in momentum space, we 
get 



H MF = E ^l D ( k )^k + N s [fi + fiK- 3J(B 2 - A' 



where 



and the matrix 



D(k) 



V k = (b kh bl H ) T 



UAlm^k) JSRe(e fc )-// 



(A12) 

(A13) 
(A14) 



where was defined in Eq. (2.5). 



The mean field Hamiltonian (A12) can be diagonalized using Bogoliubov transformation 



Hi,: "7,- \ I &fct 

-//•/■ '//, / \ b ] 



(A15) 
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After diagonalizing the Hamiltonian (A14) becomes 



Hmf = u (jlt^T + TfciTfci + l) + N s [n + fMK - 3J(B 2 - A 2 )] (A16) 
k 

The spectrum of the quasi-particle 7fc CT is 

u(k) = y/[JBRe{£ k )-n]*-[JAlm{£ k )]* (A17) 

The coefficients in Bogoliubov transformation u k and v k , which diagonalize the Hamilto- 
nian, are 

B Re(^ fc ) - H + uj k 

Z mSl (A18) 

^2(B Re(^ fc ) -fi + uj k )uj k 

Then it is easy to see that 

u k = u_ k , v k = -v_ k (A19) 

It can be shown that the above spectrum u(k) has two minimums at k — ±kc, where 
k c = (§7r,§7r). 

To describe the low energy excitation of spinon degree of freedom, the quasi-particle 7fc CT 
can be used. Consequently we are interested in its PSG. Usually the symmetry transforms 
are described in coordinate space rather than momentum space. Therefore consider the 
Fourier transformation of 7 operator at a certain momentum 

lkAr) = J2^ +q ,,e^ r (A20) 
q 

Since the low energy excitations are near ±fec, two new field operators are defined as 

q (A21) 

72a (r) = Yl l~k c +qae tq ' r 



Plug in the definition of 7^ (A15), one can get 



iqr 



Jkc+qat 

q 

Y ( U kc+q b kc+q,° + itrUko+qb-ka-9-cr) ^ ( A22 ) 
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Similarly 

72<xO) = ( y u ~kc+Q h -kc+i,<r + i°v^ kc+q b\ c __ q _)j e iqr (A23) 
q 

Another way of describing the low energy excitations of spinon field is to use two complex 
variables Z\ and z 2 , which can also be used as order parameters in the spin ordered phase. 
To begin with, z fields are introduced by discussing the ordered phase. The excitation 



spectrum (A17) has two minimums at ±&o When k becomes larger and larger the gap at 
these two momenta vanishes. Consequently the quasi-particle modes at these momenta will 
condense and form long range antiferromagnetic order. 



In the spin ordered state, spinon field condensates in the following waj^ 

izi — iz 2 \ I e lkc ' Vl 
z 2 z \ ) \e~ lkc ' ri 




(A24) 



where Z\ and z 2 are the same fields as C\ and c 2 in Ref. 031 In antiferromagnetic ordered 
phase, they acquires non-zero expectation values. However, in spin liquid phase they repre- 
sent gapped excitations. 



From Eq. (A24) one can solve z%2 in terms of spinon operators 



bit + ib], 

21 , (A25) 
bit — ib], ., 

7o _ H ikc-r 

2 -2i 

Actually the 7 operators and z fields are correlated to each other: they become the 
same thing when the system is close to the critical point between spin liquid phase and 
antiferromagnetic long range order phase. To see this, consider only the long wavelength 



component in Eq. (A21). When \q\ <C 1, we can take the approximation Uk c + q = Uk, 
Vk c +q = Vk- Therefore 71 becomes 



iq-r 



iq-T 



q 

U k C h k C +q,^ lq ' r + WV k c b ~k c -q,-a etq ' r 

q q 



Note, that 



> r ') e -i{kc+qyr' e i q -r = . ( r ) e ~ikc-r 
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and similarly 



-kc—q,—<r 



b-JrVe 



The 71 operator can be simplified as 

7i ff ~ (u kc b a (r) + iav kc bl a (r)) 



Similarly for 7 2o 



l2a ^ (u kc b a {r) - iav kc b ] _ a (r)) e 



O iko-r 



(A26) 



(A27) 



(We used Eq. flATo) ).) 

When the system is approaching the critical point (or A and B are much larger than the 



gap), from Eq. (A18) one can see that 



v k Ju kc -> 1 



(A28) 



Therefore by comparing Eq. ( A26 ), ( A27) and Eq. ( A25 ), the following relations are obtained 



7i? ~ lz i 

111 ~ ^ 

72t ~ —iz2 

721 ~ 21 



(A29) 



APPENDIX B: INTERACTION MEDIATED BY SPINON EXCITATION 

In this section we calculated the interaction between holon or doublon through exchanging 
two spinon excitations described in section |V| 

First, consider the interaction between two holons. The interaction is shown in Fig. |HJ 
in which the interaction vertex between holons and spinons comes from the hopping term 



in equation (2.3). When the spinon gap is small comparing to the band width of spinon, 



the low energy spinon excitation can be described by the Z\ and z 2 complex fields defined 
in equation ( A25[ ). In terms of z field, the interaction is 



#i„t = -tj^fifl W^i (j> Jfcc ' (rj - ri) + z* 2 {i)z 2 {j)e- ikc<rj - ri) +h. c.) +h. c. 

(ij) 



in which the terms involve high momentum transfers like snored. After 
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a Fourier transform one get from the above equation 



Him — — 2t 2J bk 2 +k z -k 1 -k A fk 1 fl 2 z{{kz)zi(k i ) Re(j^ fc4 -fc 2 + J 2 ^fc 3 +fc 2 ) 

fk 1 ft 2 z*2(k 3 )z 2 (k 4 ) Re(j 2 ^ fc4 _ fc2 + j^ 3+fc2 ) (Bl) 



k\k 2 k 3 k4 



where j is a complex number j = e 1 ^ . 



Here we are interested in the situation that both the spinon gap and the doping are small. 
When the spinon gap is small, the interaction is dominated by the spinon excitations near 
the bottom of the minimum of spinon band. When the doping is small, the Fermi surface is 
also small, so that the momentum of fermions in the interaction vertex is also small. When 
all the momenta in the interaction vertex are small, the £ function in the vertex can be 
expanded as 

& = 3 - l -k 2 a 2 + ■■■ (B2) 

where a is the lattice constant. The leading order in the expansion is a constant, which 
corresponds to a contact interaction. 

When the gap is small, the two z fields describes bosonic spinon excitations in the vicini- 
ties of momenta ± f^, The two branches of excitations has the same dispersion relation 

uj 2 k ~ m\ + c 2 k 2 (B3) 

where m s is the spinon gap and c is the spin wave velocity. The propagator of Zi fields can 



be worked out from the spinon mean field Hamiltonian (A12) through the relation between 
the two representations. The result shows bosonic excitations with the above dispersion 
relation and a spectrum weight determined by the order parameters of the spinon ansatz. 



G ^ ) = . i + ,n i + cV (B4 > 



where the spectrum weight 



z = -fi + JBxRBtkc + JAjlm^ 



Similarly the Green's function of z 2 field can be obtained. The result is the same as Eq. (B4) 



Using the interaction vertex (Bl) and the Green's function (B4) for Z1.2 fields, the Feyn 



man diagram Fig. [8jcan be evaluated. To the leading order in Eq. (B2) 

d 2 k' Z Z 



r = -*(<E/ 



(27r) 2 uj\ + m 2 s + c 2 k' 2 oo 2 + m? s + c 2 (k' — q 



2 



(B6) 
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The integral and frequency summation can be worked out at T = as in Ref. [33] 



T = -72t 2 Z 2 ^- tan" 1 ^- (B7) 
Anc 2 q 2m s v ; 

When the momenta in the interaction are all small, this is the result in the leading order. 



However, if the constraint of the complex fields in Eq. (B9) is taken into account, the contact 
interaction should give no contribution to the holon scattering process because the spinon 
fields appears in the interaction can be simply replaced by the constant 1. However, in 
the above equation we do get a contribution to the scattering amplitude from that contact 
interaction. The reason is that we use a large N expansion (or mean field approximation) 
in determining the low energy excitation of the spinon field. However, in such an expansion 



the constraint (B9) is replaced by a global chemical potential. Therefore to correctly include 
the effect of the constraint we have to take into account the 1/N order fluctuation of the 
chemical potential field. 

The original action of the z fields is a non-linear sigma model. 

£ = ^|<W (B8) 

with the constraint 

\ Zl \ 2 + \z 2 \ 2 = 1 (B9) 
The partition function of this action is 

Z = J Vz*Vz5(\z a \ 2 - 1) exp j rf D x|^ CT | 2 | (BIO) 

where in our case D — 2 + 1. In the above action, the delta function that enforces the 



constraint (B9) can be replaced by a Lagrangian multiplier field A as follows 



Z = J Vz*VzV\e^^-^ J d D x[\d^z a \ 2 + i\{\z a \ 2 - 1)]J (Bll) 

Integrate out the z a field and take the large N limit, one can get the saddle point of A field 
as iA = m 2 , where the mass is determined by the following mean field self-consistent Eq. (see 
Eq. (5.21) in Ref. 1331 ) 

rA d D k 1 1 

(B12) 



(2tt) d k 2 + m 2 g 

To calculate the correction in the 1/N order, the fluctuation of A field needs to be included. 
Section 7.2 of Ref calculated that the propagator of A field is l/EE(g), where n is the bare 
density-density correlation function of z fields: 

U(q) = [ f fc - I = ^^tan- 1 ( ^- ) (B13) 

J (27r) D k 2 + m 2 (k + q) 2 + m 2 An^fq 2 
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If the fluctuation of A field is included, there will be another Feynman diagram in addition 
to Fig. [8] (see Fig. [fijf. 




FIG. 10: Interaction 



The contribution from this diagram is 



r = (-36t 2 )n(g) • (-Uiqy^Uiq) = 3Qt 2 U(q) 



(B14) 



Therefore it cancels the contribution from Fig. 8 This is consistent with the constraint (B9) 



and it implies that there is no holon-holon scattering coming from the contact interaction. 
Consequently we need to consider interaction with momentum dependence. Taking the next 
order in interaction (|B1|), one will get 



H mt = -2t ^ $k 2 +k 3 -k 1 -k 4 fk 1 fi 2 [zi{k 3 )z 1 (k 4 ) + z* 2 (k z )z 2 {k A )] 

6-^(fc 4 -fc 2 )V-^(fc 3 + fc 2 )V 



(B15) 



where the a order gives no contribution to holon-holon interaction because of the constraint, 
and therefore can be dropped. 



Plug in the interaction in Eq. (B14), the holon-holon interaction shown in Fig. 8 and 



Fig. 10 becomes 



UAt 2 J 1 ^ 



(2vr) 2 {ul + k' 2 + niDlul + (k' - q) 2 + m 2 



d 2 k' l(k x - k'fa 2 + Kfc! + k' - qfa 2 



(2tt) 2 (oj 2 + k' 2 + m 2 s )[iol + (k' - q) 2 + m 



(B16) 
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Let ko = ki — q/2, and replace k' by k' + q/2, Eq. (B16) can be simplified as 



r = -i44tV-V 
f 3 ^ 



d 2 k' 



(k' 2 + kl) 



+ 



U4t 2 a 



2^4 



n(«) \p 



R ' 



(2tt) 2 [u; 2 + (k> + q/2) 2 + m 2 ] [uj 2 n + (k> - q/2) 2 + m 2 s ] 

d 2 k' k' 2 + k 2 

(2tt) 2 [uo 2 n + (fc' + <?/2) 2 + m 2]K + (k> - q/2f + m 2 ] 



(B17) 



It is easy to prove that the terms involving k are canceled. Therefore the interaction 
coefficient can be written as 



144tV 



R L^i 







d 2 k' 



k 



+ 



144t 2 a 4 



U(q) \/3 



-y 



(2vr) 2 [u 2 + (k' + q/2) 2 + m 2 ] [u 2 n + (k> - q/2) 2 + m 2 ] 
d 2 k' k' 2 

(2tt) 2 [uo 2 n + (fe' + g/2) 2 + ml][uol + (fe' - q/2f + m 2 ] 



(B18) 



The fc' 4 factor in the integral introduces an UV divergent. The most divergent term in the 
integral is proportional to A 3 , where the cut-off A = a -1 . However, the most divergent term 
does not depend on q and m and will be canceled out when including the contribution from 
switching the two outgoing fermions. Therefore the relevant leading order is A 2 . Combine 
that with the prefactor a 4 = A -4 , we conclude that 



A 2 



(B19) 



The dimensionality of V is [M _1 ] (see equation (B7). Therefore V should be proportional 
to the larger one of m and q when one is much larger than the other. Consequently by 
dimensional analysis we can get the results summarized in Eq. (5.5) and (5.6). 

Then let us consider the counterpart of this interaction for doublon excitations. This 
can be done in the same way as the holon interactions, while the only difference is that 
low-energy doublon excitations are around non-zero momenta fcf : (vr,0), (0, 7r) and (vr,7r). 

We start with the second term in the hopping Hamiltonian ( 2.3[ ). In the limit of small 
doping and small spinon gap, we can again replace the b ia operators by z\p fields. After a 
Fourier transform the interaction vertex becomes 



H int =2t s k 2 +k :i -k 1 -k 4 9k 1 gl 2 z* 1 (k3)z 1 (k 4 )Re(j^ k2 - k4 +j 2 ^ k2+k . i ) 

kik 2 k 3 k4 

+2t Y 5 k 2 +k 3 -k 1 -k 4 gk 1 gl 2 4( k 3)z2(k 4 ) Re(j 2 £fc 2 _ fc4 + j£k 2 +k 3 ) 

kik 2 k3ki 



(B20) 



Low energy doublon excitations are located in the vicinities of the minimums in the 
doublon band kj^x- To calculate the interaction between doublons, we need to replace the 
g operator in the above Hamiltonian with the operator g\ defined in Eq. (2.19). Then 
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FIG. 11: Interaction between spinon and holon due to the hopping term. 



Eq. (B20) becomes 



#int — 2t $k2+k3-ki-k49\k i g{k 2 Z l(k3)Zi(k4) Re(j^ fc _ DA+ fe 2 _fc 4 + j 2 ^k DX +k 2 +kz, 

kik 2 k-jki,\ 

+2^ 2J ^fe2+fc3-fci-fc 4 fl , Afcifi , AA:2 2; 2 (^3)^2(^4) R-e(j 2 ^fc DA+fc2 _ fe4 + Hk DX +k 2 +k 3 ) 
kik 2 k 3 k^ 



(B21) 



Note, that we only included interaction between the same type of doublon, because in- 
teraction between two different types of doublon cannot be mediated by spinons near the 
minimum of spinon band, and therefore is much smaller. 

With this interaction vertex and the same Green's function of z\ 2 fields we got before in 



Eq. (B4), we can evaluate the Feynman diagram Fig. 8 with incoming and outgoing type- A 



doublons together with 1/N correction in Fig. 10, The result has the same form as we have 



in Eq. (5.5) and (5.6) because the interaction vertex (B21) has the same form as the one for 
holons pi). 



APPENDIX C: BOUND STATE OF SPINON AND HOLON 



One spinon and one holon may form a bound state if we take into account the interaction 



between the two due to the hopping term in Eq. (2.3). Again, we need to include the effect 



of constraint (B9); in other words, we need to consider the contribution from the fluctuation 



of A field shown in Fig. 10 



In general, we consider the interaction between two holons and two spinons with incoming 



and outgoing momenta shown in Fig. 11 



mentum dependence according to the expansion (B2) 



The tree level interaction coming from the hopping term in equation (Bl) has the mo- 

(CI) 



H = t 



k\k 2 k-j,k4, 



fkJlX(k 3 )zAki)iG ~ a\k r + k 2 ) 2 - a\k x - fc 4 ) 



Therefore the interaction vertex at tree level is 



^(0) 
fb 



-t(6 - a 2 (fci + k 2 ) 2 - a 2 (fci - fe 4 ) 



(C2) 
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For the same reason we gave in Appendix [Bj the one loop correction due to the fluctuation 
of A field needs to be combined with the zero-th order contribution from equation CI This 
is represented by Fig. 12 and is expressed by the following equation 




v 



FIG. 12: One loop correction of the interaction in Fig. 11 



F/6 (fcl,g,*Wn) 



d 3 k' 



6-a 2 (k' 1 + k 2 ) 2 -a 2 (k' 1 -k 4 y 



1%) J (2tt)3 [„l + (k[) 2 + m ^ m + Un) 2 + (fc ; _ q) 2 + m 



(C3) 



where the interaction depends on the incoming momentum ki, momentum transfer q and 
energy transfer u n . 



Add the two contributions in Eq. (C2) and (C3) together, one immediately see that the 



terms not depend on ki and fc 3 in equation ( C2 ) will be canceled by corresponding terms in 



the other term, due to the definition of 11(g) in equation (B13). This yields to the following 
result 



r fb (kx,q,iuj n 
= 2ta 2 ki ■ fc 3 



ta 2 



d 3 k' 



2k[ ■ (k[ - q) 



(2vr) 3 K 



{k' l f + ml][{u m + u Jn Y + {k' l -q) 



mi 



(C4) 



It is easy to see that in Eq. (C4), the interaction can be seperated into two parts: the 



first term is a simple one but with a dependence on ki, which does not appear in usual 
two-body interactions; the second one is the result of a loop integral and only depends on 
the transfered momentum and energy. 

The asymptotic behaviors of the second term in m -C q and m ^> q can be studied by the 
same dimension analysis method used in Appendix [Bj By counting the power of k' in the 
loop integral one can see that the most divergent term is proportional to A, and the term 



in the next order does not depend on A. However, unlike the case in equation (B18), here 



the divergent term does give a q and u dependent interaction because of the 1 /II prefactor. 



Therefore the dominant term in equation (C4) is 



ta 



fb 



tac 2 \Jc 2 q 2 — to 2 



8nU(q,u) 



2 tan 



_^ yj c 2 q 2 —oj 2 
2m s 



(C5) 
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And it only depends on q and u. The shape of the function in Eq. (C5) is plotted in Fi 




yj c 2 q 2 — u> 2 /2m s 



FIG. 13: Shape of the interaction function in equation (C5). 



The bound state problem can be solved by solving Bethe-Salpeter equation 32 . In order 



to see the qualitative results of the binding potential in Eq. (C5) we ignored the frequency 



dependence and assume that the dispersion relations of holon and spinon are both non- 
relativistic. In this case the Bethe-Salpeter equation for the relative motion is 



jf 



— Ei 



d 2 p' 



r(p, P \K)<p( P ') 



(C6) 



2(i 7 rvr/ J (2tt) 2 

where p is the relative momentum; \i is the reduced mass in center of mass frame; E b is the 
bound state energy; K is the total momentum in the center of mass frame; <p(p) is the wave 
function of relative motion in the momentum space. Numerical calculation shows that there 



can be one or more bound states due to the interaction (C5) 



In general, the interaction depends on both the relative momentum p, p' and the total 
momentum K. However, the dominant term (C5) only depends on q = \p — p'. Moreover, 



although the full interaction (C4) depends on all momenta, the dependence on K can be 
separated 

T fb = 2ta 2 K 2 + T' fb (p,p') (C7) 

Therefore the motion of the center of mass can still be separated from the relative motion, 
and the bound state energy can be written as 



E(K) 



K 2 



2(m s + m h ) 



+ 2ta 2 K 2 + E b 



(C8) 



where E b is the bound state energy which does not depend on K, and m h , m s are masses 
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of holon and spinon respectively. Consequently the bound states have minimal energy at 
K = 0. Since the momenta we use in this Appendix are all measured from the center of 
spinon and holon energy minimums, the position of the bound state energy minimum is the 
sum of holon and spinon minimums. The former is at the origin; the later is at ± (^, f). 
Therefore the position of bound state minimums is also ± This is shown in Fig. [2] 

For doublons, we can also get an attractive interaction between them and spinon from 
Eq. ( B2 1 ) . Because of the same dimension analysis argument the interaction has the same 
form as (C5). Therefore, the bound state of doublon and spinon, or electron excitation 
locates at the sum of momenta of doublon and spinon, which is also shown in Fig. [2] 
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